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Abstract The Tutte polynomial is a powerfuU analytic tool to study the 
structure of planar graphs. In this paper, we establish some relations be- 
tween the number of clusters per bond for planar graph and its dual : these 
relations bring into play the coordination number of the graphs. The fac- 
torial moment measure of the number of clusters per bond are given using 
the derivative of the Tutte polynomial. Examples arc presented for simple 
planar graph. The cases of square, triangular, honeycomb, Archimedean 
and Laves lattices are discussed. 



1 Introduction 

It is well known that the Tutte polynomial is a powerfuU tool in order to 
explore phase transition via the Lee and Yang theorem. It allows the exact 
calculation of the partition function using Fortuin-Kasteleyn representation 
of the Potts model for a large family of graphs. This representation shows 
a strong link between phase transition and percolation purposes. We recall 
the famous paper |12j on exact critical values for site and bond percolation 
using matching and dual graph and the number of cluster per site or per 
bond. A result of Kesten is proving that bond percolation thresholds of a 
dual pair of periodic planar lattices sum 1 and then that the exact bond 
percolation threshold of the self dual square lattice is 0.5 (H]. Wierman 
[TS| establishes exact bond percolation threshold for the triangular and 
honeycomb lattice and |14) makes recent progress on bond percolation on 
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Archimedean and Laves lattices. As the chromatic polynomial is a particu- 
lar case of the Tutte polynomial, since the proof of the four colour theorem 
for planar graph a large number of conjectures have been proposed 

on the real, integer and complex zeros of the chromatic polynomial for dif- 
ferent class of graphs, see [Zj. In fact, Sokal ^l] prove that complex zeros 
of chromatic and Tutte polynomial are dense in the complex plane. A large 
class of graphs has been studied : for graphs of bounded degree we can see 
the work of Procacci et al. |1L)| and for graphs and matroids determined by 
their Tutte polynomials see and references therein. 

The paper is organized as follows. First, we recall the definition of Tutte 
polynomial of a graph and its dual. In a second part, some relations be- 
tween the number of clusters per bond for planar graph and its dual are 
established using the coordination number. The factorial moment measure 
of the number of clusters per bond are given using the derivative of the 
Tutte polynomial. The cases of self-dual graphs, square, triangular, hon- 
eycomb, Archimedean and Laves lattices are discussed in the end of this 
paper. 



2 Tutte Polynomial 

2.1 Definition 

Let be G = {V, E) an undirected planar graph, where V is the set of vertices 
and E is the set of edges. The main theorem of this paper is based on the 
Tutte polynomial of a graph G. Given a connected graph G, the Tutte 
polynomial denoted by T(G, x, y) can be obtained by computing the four 
following rules : 

(1) If G has no edges, then T{G,x,y) = 1. 

(2) If e is an edge of G that is neither a loop nor an isthmus, then 

where Gg is the graph G with the edge e deleted and G" is the graph G 
with the edge e contracted. 

(3) If e is an isthmus, then T(G, x, y) = xT{G'^,x, y) 

(4) If e is a loop, then T(G, x, y) ^ yT{G'^, x, y). 

Example : The Tutte polynomial associated with a cycle of length n, G„, 
is given by 

n-l 

T{C,,,x,y) = Y.'''+y- 
fc=i 
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2.2 Dual graph 

Given a connected planar graph G with n vertices and m edges, we assume 
that this graph has / faces. The Eulcr relation gives the following equality : 

/ = m — n + a 

with a is Euler-Poincare caracteristic. In this paper, we are interested in 
planar lattice, which corresponds to a = 1. We can define the dual graph 
G* of G : for this, it's necessary to add an exterior virtual face to the graph 
G. Thus, each vertice of G* corresponds to a face of G and two vertices 
of G* are connected by an edge if the corresponding faces in G have a 
boundary edge in common in G. We can write n* = / + 1, m* = m where 
71* is the number of vertices of G* and m* the number of edges of G*. 

Then, we have the following relation between the two Tutte polynomial 
T(G, X, y) = T{G*, y, x). For example, for self-dual graphs that means G = 
G*, the Tutte polynomial is symmetric in x and y. But, we don't have 
reciprocal result. In fact, for example we can find different planar graphs 
having the same Tutte polynomial (see Figure 1 and for details). 




Fig. 1 two graphs with the same Tutte polynomial. 

An other representation of the Tutte polynomial 

We can introduce in the Tutte polynomial the notion of rank function r : 
the rank function of a graph G is a function on the subsets of E, which 
associates to each Ac E the value r{A) = \V\ — k{G,A) where V denotes 
the set of vertices and fc(G, A) is the number of connected components of 
the graph with vertex set V and edge set A C E. This definition provides 
an explicit formula for the Tutte polynomial : 

Ace 
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In this paper, we just consider connected graph : in this case, if we set 
n = \V\, m = \E\, the Tuttc polynomial may be written as 

Ace 

^ ]k{G,A) 



(a; - 1) 

Two important values of the Tutte polynomial 

(1) If we calculate this Tutte polynomial on the curve (x — l)(y — 1) = 1 
or more specially for x = 1/p and y = 1/(1 — p), we obtain 



- "> - (.-1)1-1).. i^".'" - - fa-i).';.-i ) 

which gives 



r(G, 1/p, 1/(1 - p) = - — ^ — - — -. (1) 

(2) Now, we look at the curve {x — 1)(?/ — 1) = where q is a real. For 
X ~ 1 + 9(1 — p)/p and y — 1/(1 —p), the Tutte polynomial can be written 
as : 



Ace 



r(G, 1/p, 1/(1 -p)) ^^^^,(G.p)) (2) 



where E{q'') denote the expectation of the random variable q'' following 
the binomial distribution. This relation is well known because it gives a 
connexion between the Tutte polynomial and the partition function of the 
Potts model using the Fortuin Kasteleyn representation. 

By using this calcul, we give now a new result between the number of 
connected component for a graph G and its dual. 



3 Number of cluster per bond 

3.1 Definition and notations 

We associate with each edge of the graph G a distribution : an edge will 
be white with probability p and black with probability 1 — p. We denote 
by fc(G, A,p) the number of connected components with white edges in A, 
subset of E. 
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For the dual graph, we can build a relation between A and a set A*, 
subset of E*, the set of the edges of G* : let be e an edge of A. These edge 
is common on two faces of the graph G which correspond two vertices of 
G* and these two vertices define an edge e* in G*. So, given a subset A of 
edges e of E, the set A* of E* is the complementary set of edges e* thus 
defined. In opposite, the edge in G* shall be white with probabihty I — p 
and black with probability p. Also, if the edge e in G is white, the edge e* 
will be black in G*. We denote by k{G* , A* ,1—p) the number of connected 
components with white edges in A*, subset of E* (see Figure 2). 



Fig. 2 Example of set A in G and corresponding A* in G* . 



3.2 Connected components 

We use the following theorem and introduce a new proof based on the Tutte 
polynomial. It seems that this relation already appears in j3] in order to 
obtain some correlation duality relations for Potts model on planar graphs. 



Theorem 31 Let be G a connected planar graph and denote by G* the 
dual graph of G. 

For all subset A oi E and corresponding subset A* of E* , 

k{G*,A*, 1-p)- k{G, A,p) ^\A\-n+l (3) 
where |^| represents number of elements of A. 

In this result, we can notice that the random variable k{G* , 1 —p) — k{G,p) 
follows a binomial distribution with modality set {1 — n, . . . , n* — 1}, which 



6 



is symmetric for n = n* . 

Proof : the proof is based on the symmetry of the Tutte polynomial. For 
the graph G, on the curve (x — — 1) — q, the equation Q gives 

^ ^ ACe{G) 

and, for the dual graph, 

r(G*,l + <Z^,i)-p"'-'"*-l(l-p)l-"* y p^A-\A'\^l_p)\A'\qHG',A',l-p)-l_ 

1 — p n ^-^ 

^ A*Ce(G*) 

Notice that m* = m, to*— 7i* + l = /* = 1 and 1— 7i* = — / = —m+n—l, 
so that the previous equation can be written as 



^ ^ A*Ce(G*) 

(4) 

Moreover, the following relation T{G*,x,y) ~ T{G,y,x) gives us an other 
way to calculate (@J, namely : 

T(GM + g^,i) = r(G,i,l + 9T^) 



p_\m / „ p \ / 1 \ 



1- 
P 



1— p ^ 



It becomes 

r(G*,l+q^,i)=/-"(l-p)"-l y pl^l(l-p)-|^lgWG,A,p) + |A|-n) 
' ' ACe(G) 

(5) 

Before, we made correspondence between the set A and the set A* , namely 
an edge e {e G A) is white in G when the corresponding edge e* is black in 
G*. So, we conclude |A| + |^*| = to. A* can be seen as the complementary 
of ^ in G*. In this way, the equation gives us 

r(G*,l + <Z^,i)=pl-"(l-pr-l y pl-4|(l-p)-|^lg'=(G*,AM-rt-l^ 
^ ^ ACe(G) 

Finally, if we compare this equality and the relation it becomes 
k{G*,A*,l~p)- k{G,A,p) = \A\ -n+1. 
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3.3 Mean of connected components 

The previous result yields a quick way to compare the means of k{G* , 1 —p) 
and k{G,p) under the uniform distribution : 

Theorem 32 

E(k(G*, 1 - p)) - E(k(G, p)) = mp - n + 1 

In fact, by multiplying equality JSJ by p'"'^' (1 — p)™^''^' and by summing 
on all subsets A oi E, we get : 

^ p\^\\-pr-\^\k{G\A\i-p)- p\^\\-pr-\^\k{G,A,p) = 

ACe(G) ylCe(G) 

^ pl^l(l-pr-l^l(|A|-n+l). 

ACe(G) 

By the correspondence between the sets A and A* and the relation \A\ + 
\A*\ — m, this equality becomes : 

^ p™-i^*i(i-p)i^*ifc(G*,^M-p)- p'^'(i-pr"'^'fc(G,A,p) = 

A*Ce(G*) ACe(G) 

5]0'(i-pr-'(/-n+i). 

which gives 

E(k(G*, 1 - p)) - E(k(G,p)) = mp - n + 1. 

Some remarks : (1) From the theorem 1321 we can notice the particular 
status played by the value p = {n—l)/m. For this value, the two means are 
equal. We'll discuss on this value in the next section for different graphs. 
(2) Following [H], we can compare the expected values of the length of the 
minimal spanning tree of G denoted by E{Lmst{G)) and of G* denoted 
by E{Lmst{G*)). They prove the following relation 

E{Lmst{G))^ f E{k{G,p))dp-l. 
Jo 

In our context, we obtain 

m 

Ei-p(Lmst(G*)) - Ep(Lmst(G)) = - - n + 1. 

In this equality, it appears the value m = 2{n — 1). This relation between 
number of edges and vertices of G is a necessary condition so that the Tutte 
polynomial may be symmetric in x and y. 
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3.4 An other approach 



The result of theorem 1321 can be presented under another form : if we go 
back to the equahty (0, we infer that 



gr(G,l + g<i^,^) 

t{gMpM{i-p)) 



and 



k(G*.l~p)\ _ ^ ^ J I 'i i^pi pj _ -rt \ ; pJ ' 'i i-pj 

' " Tie 1/p, 1/(1 - p)) " T(G, 1/p, 1/(1 - p)) ■ 
These two expressions are differentiable in q so that 



i?(fc(G,p)g'=(^^^')-i) 
and 



E{k{G\l~p)q'''^^*'^~P^-^) 
Takinq 1 for q, we get : 



r(G, 1/p, 1/(1 -p)) 
r(G,i,i + q^ 



+q- 



p i-p 



T(G, 1/p, 1/(1 -p)) 
P> P dy p 1-p 



T{G,i/p,i/{i-p)) '1-P r(G, 1/p, 1/(1 -rt) ■ 



EikiG,p)yEik{G\l^p))^pil-p) 
which can be written as 



^ dx^^^' p' 1-p^ {I-pY dy^^'^' p' l-p\ 



r(G, 1/p, 1/(1 -p)) 



Proposition 31 

E{k{G,p)) - E{k{G\ 1 -p)) = -p(l - p)^\ogT{G, -^). (6) 

dp p 1 — p 

By using the value Q), we recover the relation of the theorem. 
Remarks : 



1) For self dual graphs (G = G*), it is natural to take p = 1/2 in 
©. For this value the right term in © vanishes, thus p = 1/2 can be 
seen as a minimum of the Tutte polynomial T{G, x, y) on the hyperbola 
{x — l){y — 1) = g = 1 with x = 1/p and y = 1/(1 — p). We can sec that 
p — 1/2 = ■y/g/(l + ^/q) for q = 1. For other real values of q, we have a 
similar result that p = y^/ (1 + ^) is a minimum of the Tutte polynomial 



9 



on the hyperbola {x — l)(y — 1) = q. Notice that this particular value is 
conjectured to be the critical probability for q-Potts model. 

In fact the Tutte polynomial for self dual graphs is symmetric and can 
be written as 

n — 1 a 

where n is the number of vertices of G and Oq.^ are the coefficients of 
T{G, X, y). On the hyperbola (x = 1 + q/v^ y = \ + v), 

— = -q/v^ dy (foi^ ^ 3 (Py ^ ^ 

ris — l , 

We can always write T{G, x,y) ^ V cJx+y)°' = T(G, v) thus —T(G, v) 

^-^ dv 

a=0 

can be factorized by 1 — q/v'^ and vanishes for = q which corresponds 
to X = y. By showing that the Tutte polynomial is a convex function of v, 
V > otherwise p g]0, 1[ , v = ,Jq reaches the minimum and is equal to 
T(G,1 + + 

2) From the equation we obtain 

^ t{gMpM{\-p))- 

By i differentiations with respect to and taking g = 1, the factorial 
moments are 

. ^\t{c\ ^ ) 
V» e {1, . . . , n}, E\{k - l)(fc - 2) . . . (fc - »)] ^ (1^)' >' ; 

(7) 

9^ . 

where tt—T is the ith partial derivative with respect to x of Tutte poly- 
ox* 
nomial. 



4 Applications 

4.1 Finite planar graphs 

As example of finite planar graphs, we study the wheel which is a self 
dual graph. The wheel is the graph obtained from a cycle by adding a new 
vertex adjacent to all vertices in the cycle. This graph contains n vertices 
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and 2(71—1) edges. In this case, the equation in theoremK^2lcan be written 
as follows 

E [k{G\ l-p)]-E [k{G,p)] = 2{n - l){p - i) (8) 

The natural choice of p is 0.5 and so that for all graphs with m = 
2(n — 1). This condition is necessary for having a symmetric Tutte polyno- 
mial and in particular this is satisfied for self dual graphs. 

Now, by looking at the result of the theorem 1321 and the previous ex- 
ample, we are led to study the particular value p^ = -^^^ for general finite 
planar graphs. It comes 

n 

Y,E{\C,\-') - E{k{G,Pe)) = E{k{G*, 1-pe)) 

where Gi is the cluster containing the site i. 

In table we present some graphs with their particular value p^. For 
two firsts graphs, s denotes respectively the number of squares and trian- 
gles. The third graph is the square lattice at step s = 2. For the three 
following graphs, we just give the pattern associated with the square lat- 
tice, respectively : bow-tie lattice, crossed matching and octagonal. The two 
last graphs are respectively the honeycomb lattice and the kagome lattice 
at step s = 2. For all graphs, we compute, at each step s, the number of 
edges and vertices and the corresponding value Pe- 

Remark that 

_ n - 1 _ 2(n - 1) 
where qi denotes the number of neighbours of the site i. 



4.2 Infinite planar graphs 

Taking the thermodynamic limit in pe, we obtain pi = ^ where q is the 
coordination number of the graph. Obviously from Euler relation 

^ n - 1 n* — 1 

Pe+Pe ^ \ = 1. 

m m 

Of course, for self dual graph, pi = 0.5 : it corresponds to the critical bond 
percolation threshold of the square lattice. 
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Name 


Graph 




1^1 


V -i 

- \E\ 


lim pe 

s — > + oo 


s squares 




O „ 1 1 

OS + i 


2s + 2 


2s + 1 
3s + 1 


Zj 6 














s triangles 




2s + 1 


s + 2 


s + 1 
2s + 1 


1 /o 

1/2 






Square Lattice 




2 


2 


2 




4s(2s - 1) 


(2s)2 


2s + 1 
4s 


1/2 




2 


1 


2 








2 


2 






Bow-tie Lattice 


C2 = 








(4s - 2) (10s - 3) 


(4s - If 


4s 


2/5 


\ 




10s -3 


Crossed Matching 


C3 = 








4(4s - 2) (3s - 1) 


(4s - 1)^ 


s 


1/3 






3s - 1 


Octagonal 








4(2s - l)(3s - 1) 


8s^ - 4s + 1 


s 

3s - 1 


1/3 


Honeycomb Lattice 


\ ^ ) — \ ' ) — \ ^ / 

\ / 2 \ ( 2 \ / 


12s^ - 4s - 2 


8s^ - 2 


Ss^ -3 
12s2 - 4s - 2 


2/3 


kagome lattice 




2(2s - l)(6s - 1) 


4s(3s - 1) 


(6s + 1) 
2(6s - 1) 


1/2 



Table 1 Some graphs with their particular value Pe- 
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111 order to understand the last column of the table 14. l1 we write down 
the equation of theorem 1321 as follows 

E(k(G*, 1 - Pe)) - E(k(G, Pe)) = mpe - n + 1. 

By taking pe — ^7^1 we obtain 

E(k(G*, 1 - ^)) = E(k(G, ^)). 

Ill 111 

When n tends to infinity with m, m function of n. it becomes 
E(k(G*,l- lini pe))=E(k(G, lim pe). 

n — * + oo n — »+oc 

The last column of the table l4.ll gives value of these limits for our particular 
graphs. For honeycomb lattice G and its dual, the triangular lattice, we find 

E{k{G,2/3)) ^ E{k{G*,l/3)). 

How to understand for example this 2/3 associated with honeycomb lat- 
tice ? We can break up the expression of Pe as 

_ 2(n-l) _ 2 

E * + E ^' 7 E + 7 E 

— ' — ' n — 1 ^ — ' n — 1 ^ — ' 

with Hint + Tiext = n. Ugxt IS the number of vertices of the convex hull 
associated with honeycomb lattice at step s. In this case and for all regular 

lattice with bounded degree, lim„^_(_oo — = 0, which gives 

n 

y y 2(n-l) 2 
lim Pe = hm — — ^ -. 

i=l 

So that to calculate this limit, we avoid boundary effects. We only need of 
the coordination number of the graph : 3 for the honeycomb lattice. So we 
can state for Archimedean and Laves lattices the following relation 

E{k{G,-_))^E{k{G\l --_)). 

Concluding remarks 

1) The arithmetic and harmonic means of the number of neighbours of 
each vertex of the graph are the same for Archimedean lattices and 
also for Laves lattices with different weights as explained in Q]. In this 
paper, he builts a new mean with other weights in order to classify 
graphs with different bond percolation thresholds. Otherwise, results 
about ordering critical site and bond values with different connective 
constant are presented for some Archimedean lattices in |14j . 
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2) 111 the preceding decomposition of q, one may ask what is the role of 
part concerning the vertices of the convex hull, it remembers of course 
conjecture connecting isoperimetric constant to percolation threshold 
of a graph. For example, it could be instructive to compare the speed 
of convergence of Pe to | for different graphs. 

3) Chang and Shrock ^ remark that for strips of graphs such as square, 
triangular, kagome, that complex singularity of the number of cluters 
per bond are precisely on the boundary of the accumulation set of the 
zeros of the partition function. 

4) The present study shows the major role of self-dual graphs. An inter- 
esting perspective will be the location of complex zeros of their Tutte 
polynomial and a discussion on the thermodynamic questions. 
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